Distribution of NMR relaxations in a random Heisenberg chain 
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Nuclear magnetic resonance (NMR) measurements of the 29 Si spin-lattice relaxation time Ti were 
used to probe the spin-!/2 random Heisenberg chain compound BaCu2(Sii- a: Ge I )207. Remarkable 
differences between the pure (x = 0) and the fully random [x = 0.5) case are observed, indicating 
that randomness generates a distribution of local magnetic relaxations. This distribution, which 
is reflected in a stretched exponential NMR relaxation, exhibits a progressive broadening with de- 
creasing temperature, caused by a growing inequivalence of magnetic sites. Compelling independent 
evidence for the influence of randomness is also obtained from magnetization data and Monte Carlo 
calculations. These results suggest the formation of random-singlet states in this class of materials, 
as previously predicted by theory. 

PACS numbers: 75.50. Lk, 76.60.-k, 75.40.Gb, 75.10.Pq 
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Random variations of the exchange coupling constant 
J in ferro- and antiferromagnetic systems (bond random- 
ness) can have a profound effect on their magnetic prop- 
erties. An innovative scheme for studying random spin- 
x ji Heisenberg chains (RHCs) was developed already in 
1979 by Dasgupta and Ma 0. Their real-space renormal- 
ization group (RSRG) method was extended by Fisher 
Q , and later employed for a much larger variety of prob- 
lems related to quenched disorder in quantum magnets. 
Among these are chains with random ferro- and antifer- 
romagnetic couplings disordered spin-1 chains and 
spin-!/2 ladders 0Jp, Heisenberg magnets in two and 
three dimensions |6[, as well as dilute spins in doped 
semiconductors Q. At the core of the RSRG method is 
the so-called decimation procedure. For RHCs it involves 
an iterative suppression of the degrees of freedom of the 
most strongly coupled spins via the formation of "frozen" 
singlets. Fisher [2j has shown that this procedure con- 
verges to a universal fixed point. The resulting ground 
state, characterized by spins coupled at all possible dis- 
tances and energy scales, is also known as the random 
singlet (RS) state. At low temperatures, only excita- 
tions of the most weakly bound singlets are expected to 
be relevant Q. As a result, regardless of the microscopic 
details of the disorder, a large class of bond-disordered 
systems share universal low-temperature and low-energy 
properties. In recent years, besides providing predictions 
for bulk properties, new theoretical work has started to 
address also the problem of dynamical observables || [t| . 

To date only a few good realizations of RHCs have 
been achieved. Some early examples f]~oj ] were later 
found to have additional intricacies beyond the Heisen- 
berg chain model 
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Usually RHCs are obtained as 
solid solutions of two prototypical spin-chain compounds 
with different physical properties. A well known exam- 



ple is Sr 3 CuPti_. ! Tr. E 06 [12j . whose bulk properties are 
in agreement with the RSRG predictions. To obtain a 



direct microscopic insight on the problem, we employed 
NMR experiments to study the local magnetic relax- 
ations in a similar system, namely BaCu 2 (Sio.5Ge . 5)207 
[l^ . In this Letter we show that a strongly temperature- 
dependent distribution of magnetically inequivalent local 
relaxations indeed exists in an RHC material. In the con- 
text of the above-mentioned decimation procedure, we 
regard this phenomenon as an indication of the forma- 
tion of an RS state. 

The material BaCu2(Sio.sGeo. 5)207 considered here, 
is isostructural with its parent Si- and Ge- compounds. 
The quasi-lD system BaC^SiaO? (see Fig. [1]) has an 
in-chain coupling J = 24.1 mcV II 41 and orders anti- 
ferromagnetically at Tn = 9.2 K [15j. Bond random- 
ness in BaCu2(Sio. sGeo. 5)207 is introduced via the dif- 
ferent bond angles of Cu-0 exchange paths, a direct con- 
sequence of the different Si and Gc covalcnt radii. A 
previous study [Tol ] revealed a logarithmic temperature 
dependence of the spin susceptibility x(^) &t l° w temper- 
atures [l[, a characteristic RHC feature predicted to be 
independent of the initial J distribution. Despite that, 
and following some original misinterpretation, inelastic 
neutron scattering experiments were found to be quan- 
titatively consistent with Luttingcr-liquid behavior ex- 
pected for a disorder- free spin chain with a single effective 
coupling J ofr = (J Si + Jgc)/2 « 37 mcV [IB, tl3l ■ lt was 
concluded that the RS-relatcd physics manifests itself at 
energies much lower than those accessed by neutron stud- 
ies. Probing these lower-energy/longer-time scales with 
NMR was the main motivation of the present work. 

High-quality single crystals of BaCu2(Sii_. T Gc x )207 
with x = and x = 0.5 were grown by the floating zone 
technique. The magnetic susceptibility was measured via 
dc SQUID magnetometry. Spin-lattice relaxation times 
Ti were measured using standard pulse techniques. Due 
to the rather long relaxation times (tens of seconds at low 
temperatures), an aperiodic saturation recovery with an 
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echo detection was chosen. With a field (IqH = 7.024 T 
applied along the a-axis, the 29 Si NMR signal was found 
at 59.414 MHz. 
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FIG. 1. (Color online) Crystal structure of BaCu2Si207 em- 
phasizing the copper chains (dashed lines) and the NN con- 
figuration of the silicon atoms. 
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FIG. 2. (Color online) Magnetization M per Cu ion of BaCu2- 
(Sio.sGeo. 5)207 as a function of temperature (a), and field (b), 
normalized with respect to the saturation moment, assuming 
g L = 2. The fit of M(T) with an RHC model was adopted 
from Ref. i, while the M(T) and M{H) curves for chains 
with a single isotropic exchange value were calculated as in 
Ref. [Isl. QMC simulations of RHCs are also shown (see text). 
The symbols are explained in the panels. 

That it is precisely the low-energy dynamics of BaCu2- 
(Sio.5Geo.5)2C>7 which is most affected by randomness 
has a clear confirmation in the measured bulk magnetic 
response (Fig. [5]). Unlike for disorder- free spin chains 
(light dashed lines), the magnetization M(T) of the dis- 
ordered material (dots) tends to diverge at low temper- 
atures. This divergence is notably different from that 
observed for x — 1, where it is related to weak ferro- 
magnetism and thus strongly orientation-dependent fl9| . 
In our case, instead, the M(T) curves only reflect the 
weak anisotropy of the g-factor tensor. We interpret 
the x = 0.5, low-T magnetization by the RS theory, 
which attributes it to the low-energy states associated 
with the above-mentioned weakly-bound singlets. It prc- 



to the data above 4 K with a cutoff ft ~ 66.3 ± 0.7 
meV (dark dashed line). Since the iterative decimation 
procedure suppresses the strongly bound singlets with 
an effective coupling larger than T, the fraction of un- 
paired spins is tit ~ l/[ln(X2/T)] 2 . The paramagnetism 
of this fraction implies x(T) ~ tlt/T, with £1 the tem- 
perature above which the entire chain is effectively in a 
paramagnetic regime. The characteristic energy scale of 
the disorder-induced excitations is, however, much lower 
than f], because already a 5 T magnetic field seems to 
saturate them. This is clearly evidenced by the slope 
of M(H, 2K), which at higher fields tends to that of 
a disorder- free chain (Fig. [2b)- Quantum Monte Carlo 
(QMC) simulations of Heisenbcrg spin- 1 /^ chains, using 
6000 spin sites with randomly distributed but equally 
probable Jsi and Jg c couplings, confirm these results. 
M(T) and M(H) averaged over many (JV > 40) ran- 
dom configurations were obtained via a directed-loop al- 
gorithm [iol within the ALPS 2.0 package [2l| and are 
shown in Fig. [2] The M(H) comparison is made at 10 
K, where the system can be treated as ID. Considering 
the lack of any free parameters, the agreement with data 
is remarkably good. Incidentally, a fit of the RS predic- 
tion to the low-T QMC results provides il Jq c , thus 
reinforcing the pertinence of the RS theory. 

Similarly to the magnetization data, also 29 Si NMR re- 
laxation rates T^ 1 for x = and x = 0.5 show remarkable 
differences (see Fig. [3]). The local magnetic field at the sil- 
icon nuclei originates from the Cu electronic spins Sj and 
the relevant hyperfine interactions (see Fig. [TJ . Provided 
a matrix description of the latter is adequate and the 
hyperfine coupling tensor at the j-site for the c-chain is 
A>' c (see Fig.[T]for the index notation), the g-space quan- 



tities A c aJ3 {q) 



with a, j3 = {x,y,z} 



determine the longitudinal relaxation rate: 



±- K h 2 k B T^ [A 2 xa {q) + A 2 ya {q) 



-x,y,z 



0J„ 



(1) 



diets x(T) ~ T-Mrr^/T) [2j, providing a good fit 



where u n is the resonance frequency. Detailed calcula- 
tions reported elsewhere [22| show that the q-dependent 
prefactor in Eq. ([1}, also termed the NMR form factor. 
can be expressed as f 2 (q) = M a + P a cos [q ■ (r 2 — 
with M a and P a as constants. 

The q-sum in ([1]) was evaluated numerically using the 
dynamical susceptibility of an isotropic spin-Y2 Heisen- 
berg chain given by theory [23} . Because of the large ex- 
change coupling (J = 24.1 meV ^> /xb-B ~ 0.4 meV), spin 
chains in BaCu2Si207 were considered as non-interacting 
above the ordering temperature and effectively in zero 
field. Employing the simplest linear description of a Lut- 
tinger liquid, i.e., without any scattering of quasiparticles 
and neglecting the influence of the form factor, the domi- 
nant contribution in ([T]) is at q = it and T\ is temperature 
independent. It turns out, however, that including f 2 (q) 
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FIG. 3. (Color online) 29 Si NMR spin-lattice relaxation data of BaCu2(Sii_ I Ge a; )207 for x = and x = 0.5 (no and maximum 
randomness, respectively), a) relaxation rates for both cases, b), c) recovery of magnetization, and d) T-dependence of /3 for 
x = 0.5. The symbols are explained in the panels. 



is essential to obtain a reasonable fit to the T-j - data for 
x = 0. Considering results in an increase of T^ 1 

with T, as observed and shown in Fig. [3] A similar in- 
crease, without the influence of may also be due 
to a non vanishing scattering of quasiparticles, implying 
that the dominant contribution is at q = 0, as shown in 
Ref. 12J. A convincing demonstration of the different im- 
plications of q = and q = 7r excitations on the tempera- 
ture variations of T-T was given in Ref. [25j for Sr 2 Cu03. 



Additional data of with different crystal orientations 
(not shown here because of space restrictions) confirm 
that fa(q) cannot be neglected in our case. This fact 
and the verified absence of significant spin-diffusion ef- 
fects on the field dependence of T^ 1 indicate that the 
dominant component is at q = ir. In Fig. [3] we display 
the experimental results together with fits considering 
the two possibilities mentioned above. Including 
and presuming the zero-field limit, only J, M, and P 
(the latter two independent of a), are required as free 
parameters for the fits. The best results are obtained 
with J ~ 55 meV, distinctly larger than expected. We 
recall, however, that the evaluation of J via NMR is not 
straightforward. It enters via a non-zero form factor and 
its uncertainty depends on the choice of the temperature 
range of the fit. Nevertheless, the results for the x = 
case show that 29 Si NMR clearly probes the Luttinger- 
liquid physics, hence providing a convenient technique for 
studying the low-energy excitations also for x > 0. 

Referring to the latter case, we notice that due to 
the loss of translational invariance in disordered systems, 
Eq. ([1} is not valid. The translational inequivalence of 
sites introduces a multitude of local nuclear relaxation 
times. For a nuclear spin / = !/2, the individual relax- 
ations are simply exponential, but collectively they give 
rise to a system-dependent and generally non-exponential 
relaxation. If r is one of the local relaxation times, a 
possible interpretation for the observed stretched expo- 



nential recovery is provided by 

o 

p(r,T) e _T <^> dr = e~\-~ 



~i . (2) 



This equation implies that if the same type of relaxation 
occurs with a probability p(r, T) and with a characteris- 
tic temperature-dependent r(T) at different nuclear sites, 
then the global relaxation function is a stretched expo- 
nential, characterized by the parameters /3(T) and t s (T). 
Relaxations described by a stretched exponential are 
quite common in disordered magnets. The temperature 
dependence of the exponent j3 has been reported for a 
number of different systems j^, such as spin glasses, 
quasicrystals, polymers, etc. In our case, while for x = 
a simple exponential is adequate for describing the mag- 
netization recovery, stretched exponentials are needed for 
x = 0.5. The corresponding data and fits are shown in 
Fig- Eh and|3t, respectively. From the latter, t s {T) = T\ 
and /3(T) are obtained and shown in Figs. [3^ and[3Ji. 

Extracting the weights p{r 1 T) via an inversion of 
Eq. ([2]) is a mathematically non-trivial task 27j. To ac- 
count for our case, where (3 > 0.5, we start from Eq. ([2]). 
but now using the variable s = t s /t j29j: 



<t/r. 



P(s,/3) e- st/Ts ds, 



(3) 



where P(s,/3), as p{r) before, represents a particular 
probability distribution function (PDF) of relaxation 
rates for a given value of /?. The inverse of the Laplace 
transform of ([3]) [29[ provides the required relation be- 
tween P and p, the probability distributions of relaxation 
rates and times, respectively: 



T / T 



2 ' 



(4) 



Although P(s,j3) is obtained only numerically, Eq. (|4]) 
allows us to extract the distribution of the local relax- 
ation times in a random Hcisenberg chain starting from 
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the experimental NMR data. Results of this procedure 
for a scries of different temperatures are shown in Fig. [4] 
At high temperatures, where j3 — > 1, the distribution 
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FIG. 4. (Color online) Probability distribution of NMR spin- 
lattice relaxation times in BaCu2(Sio.sGeo. 5)207 at different 
temperatures, as derived from fits to the experimental data 
partly shown in Fig. [3J:. 

of relaxation times tends to a delta function peaked at 
t = t s (T), while at low temperatures the distribution 
diverges weakly as r — > and hence remains integrable 
27| . In all cases the distribution is normalized and does 
not depend on temperature, thus reflecting the conserved 
total number of magnetic sites. In a disordered system 
the sites are magnetically inequivalent, giving rise to a 
broadening of the relaxations' distribution at low tem- 
peratures. For T — > 0, the distribution p(r, /3) peaks at 
t — ?• 0, but its mean value t s increases due to the longer 
tail of the distribution. 

We argue that the measured NMR parameters can be 
used to obtain information on the low-energy dynam- 
ics of RHC systems, inaccessible by standard neutron- 
scattering techniques [HI]. In case of randomness, the 
lack of translational invariance prompts for models based 
on real- instead of reciprocal-space configurations. In the 
RS picture, upon lowering the temperature there is an in- 
crease in the number of singlets with different couplings 
being formed, while nx goes to zero. This results in a 
continuous distribution of magnetically inequivalent sites 
and in longer t s (T) at lower temperatures, as reflected 
in our data. In this low-T scenario the low-energy dy- 
namics is determined by the increasing importance of 
the "frozen" singlets, at the expense of the yet uncou- 
pled spins. We shall address this point in a forthcoming 
publication [22| . 

In conclusion, we have carried out a detailed 29 Si 
NMR study of the spin-lattice relaxation rates in the 
BaCu2(Sii_ a: Ge a ;)207 system for x = and x — 0.5, cor- 
responding to no and to maximum disorder, respectively. 
NMR turns out to be an appropriate technique for the 
study of low-energy dynamics in gaplcss antiferromagnets 



with bond randomness. Guided by the experimental ob- 
servation that the relaxation of magnetization is charac- 
terized by a stretched exponential, we show that it is pos- 
sible to extract the distribution of relaxation functions in 
a random Heisenberg chain from experiment. With the 
additional evidence from magnetization data and QMC 
calculations, we interpret this result to reflect the for- 
mation of random-singlet states in a random Heisenberg 
chain, as predicted by theory. 
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and the NCCR research pool MaNEP of SNF. 
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